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Abstract The starburst IMF is probably as much of theoretical interest and prac- 
tical relevance as it is a subject of observational controversy. In this 
conference paper we review the most common methods (star counts, 
dynamical masses, and line ratios) to derive or constrain the IMF, and 
discuss potential problems and shortcomings that often lead to claims 
of anomalous IMF slopes and mass cut-offs. 



1. Introduction 

The IMF in regions of violent star formation activity, so-called star- 
bursts, is still a rather controversial issue. For instance, the ultra- 
luminous infrared galaxy Arp 220 with a total gas mass of 1.6 x 1O 1O M 
Arp 220 has about five times the amount of gas available in the Milky 
Way. However, its star formation rate of 300MQ/yr is about 200 times 
higher than in our Galaxy, and - in addition - the regions of star for- 
mation are rather concentrated near the center(s) of the galaxy. 

It may seem plausible that in such regions, where the physical bound- 
ary conditions for star formation differ so substantially from the solar 
neighborhood, the IMF may be different as well. Such a difference may 
be indicated by either a different IMF slope ( "top heavy" ) , a significantly 
higher low- mass cut-off ("2 — 5M ") or a different high-mass cut-off. 

The observational verification of these differences is all but trivial: 
there are few regions that can serve as starburst templates, most of 
which are at large distances, have unknown morpholgies, and are often 
heavily extincted. For closer regions, which can be resolved into individ- 
ual stars, the dynamical range of typically more than 10 magnitudes in 
flux differences between low and high mass stars is a significant limita- 
tion. Not surprisingly, the strongest support for anomalous IMFs to date 
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comes from the observations of distant starburst galaxies or individual 
super star clusters. While those environments are indeed rather different 
from the local ones it should also be kept in mind how little we know 
about these distant system as the uncertainties in both observations and 
models are fairly large. In this conference paper we discuss the general 
difficulties in the IMF determination. 

2. Method I: star counts 

The most common method, which has also be shown to work well 
in nearby regions of quiescent star formation, is the direct detection of 
stars. After photometric extraction of the sources and subtraction of 
foreground and background field stars, the correction for incompleteness 
(see below) takes care of instrumental limitations and possibly of the 
extra contribution from unresolved binary systems. Finally, the conver- 
sion from magnitudes to masses leads to a present-day mass function, 
which is often, and not very precisely, called the initial mass function as 
dynamical effects in young clusters are (incorrectly) considered to be of 
little importance. 

The results from a recent study of the IMF in R136, the ionizing 
cluster in the 30 Doradus nebula, are shown in Figure 1 as a function 
of radial distance (Andersen, Brandl & Zinnecker 2004). The errors for 
the innermost 1 pc (4") are very large and the derived slopes not very 
meaningful. However, outside the core region the slope agrees, within the 
uncertainties, well with a Salpeter (1955) slope and shows no significant 
steepening with distance to the cluster center. 

However, each step from the observed image to an IMF slope involves 
potential difficulties, of which we can here address only a few: 

2.1 Problem: the high- mass stars 

At no time within a massive, coeval cluster's lifetime will it be pos- 
sible to observe all its members on the main sequence. Many of the 
most massive stars will have collapsed or exploded before the lower mass 
stars reach the main sequence. High mass stars show strong luminos- 
ity evolution and mass loss during their short lifetimes. In particular 
the luminosity evolution of a given spectral subtype within its lifetime 
can be larger than the initial difference between different spectral sub- 
types. It is therefore almost impossible to accurately determine masses 
of O-stars from photometry alone (and individual spectra are usually 
not available). 

There is also mounting evidence that the evolution of massive stars in 
the core of a dense stellar cluster is significantly influenced by the other 
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Figure 1. The derived IMFs for R136 in radial bins using a Siess et al. (2000) 3 Myr 
isochrone with half-solar metallicity, a distance modulus of 18.5 and an extinction of 
A v = 2 (Andersen, Brandl & Zinnecker 2004). 



cluster members. Portegies Zwart et al. (1999) have shown in N-body 
simulations that due to the high stellar densities mergers of massive 
stars occur frequently, and by doing so increase its cross-section while 
the cluster is becoming denser due to mass segregation, in some kind of 
"runaway merger" process. In other words, a central massive star grows 
steadily in mass through mergers with other stars in less than 3 — 4 Myr, 
and the observed properties are not representative of the initial mass 
distributions. 

2.2 Problem: the low-mass stars 

Besides the obvious difficulty of detecting relatively faint, low-mass 
stars a significant uncertainty will be introduced by the fact that in 
most young clusters they are still in their pre-main sequence (PMS) 
phase. (The fact that PMS stars are brighter than their main sequence 
counterparts at infrared wavelengths, improves the situation slightly). 
However, the models for converting magnitudes to masses depend crit- 
ically on the cluster age which has typical uncertainties in the order of 
0.5 Myr. Usually all stars within a cluster are assumed to be co-eval. 
Furthermore, the conversion also depends on the wavebands that were 
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used and on the specific PMS model. This situation is illustrated in 
Table 1. 

Table 1. Comparison of IMF slopes derived from the same data but with different 
PMS models, age assumptions, distance moduli or based on different filter bands. 
The differences in the resulting slopes are striking! Table from Andersen (2004). 



Model a 


Age [Myr] 


DM" 


color 


slope 


SDF 


1 


13.9 


Js-Ks 


-2.20 ±0.13 


SDF 


1 


13.9 


H — K 8 


-1.46 ±0.23 


SDF 


1 


14.3 


Js-K s 


-2.41 ±0.14 


SDF 


2 


13.9 


H — K s 


-1.46 ±0.34 


PS 


I 


13.9 


Js-K s 


-1.70 ±0.19 



"SDF - Siess et al. (2000); PS - Palla & Stahlcr (1999). 'DM - distance modulus 

The effects of heating, outflows and stellar winds (Lamers, Snow 
k, Lindholm 1995), and in particular the strong winds from O-stars 
(Churchwell 1997), will cause the residual gas to be removed from the 
cluster on rather short timescales of 1 Myr. If the gas cloud was orig- 
inally virialized, and the star formation efficiency is about 30%, the 
cluster will loose about 2/3 of its mass. As a result, a signifcant fraction 
of its young members will get lost, preferentially low-mass stars at the 
high tail of the velocity dispersion and/or at large distances (Kroupa <fe 
Boily 2002). 

This has important observational consequences: even for a rather dis- 
tant Galactic cluster like NGC 3603 at 7 kpc with an assumed velocity 
dispersion of 5km/s the radial distances of initial cluster stars can be 
as large as 100" after 1 Myr (including a cos(45°) projection factor). If 
the reference field taken to estimate the field star population is closer 
than this radius there are two problems: (i) stars that were born in the 
center are now missing there, and (ii) those stars may now be mistaken 
as field stars and subtracted from the cluster sample. Both effects work 
in the same direction and introduce an artificial flattening of the IMF. 

2.3 Problem: differential extinction 

Sirianni et al. (2000) reported a "definitive flattening below « 2M Q " 
of the IMF slope at r > 1.5 ps around R136. This finding was based 
on optical data from HST/WFPC2 and the numbers were corrected for 
incompleteness. Unfortunately, the standard Monte-Carlo technique for 
incompleteness correcton - adding artificial stars to the "real" image 
and computing their detection probability as a function of distance from 
the cluster center and source brightness - corrects only for statistical 
(crowding, blending) and instrumental (noise, flatfield) effects. In con- 
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trast, strong variations in the more localized differential extinction are 
caused either by the evolution of the PMS object or by rather compact 
and patchy foreground extinction, and can vary by several magnitudes 
at optical wavelengths. These variations cannot be easily corrected with 
Monte-Carlo techniques (see Figure 2). Since the reddening distribu- 
tion is very uncertain, the best work-around are observations at near-IR 
wavelength where the extinction is reduced by about an order of magi- 
tude. 
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Figure 2. The influence of differen- 
tial extinction on the observed IMF 
slope. For a constant detection limit 
(dashed line) variable reddening can 
cause an artificial "flattening" . Figure 
from Andersen (2004). 



Figure 3. The effect of spatial reso- 
lution on the derived half-light radius 
r M , measured for R136 locally, and re- 
measured at larger distances (Brandl 
et al. 2005). 



3. Method II: cluster mass-to-light ratios 

Most starburst systems are located well beyond the Local Group and 
are not resolvable into their individual stellar members (except for the 
very brightest stars) with current techniques. However, so-called super 
star clusters (SSCs) may be the result of a recent starburst, and their 
integrated properties can be observationally studied. 

3.1 The basic principle 

For a given set of assumptions (see below) the dynamical mass M<± yn 
of a cluster is related to its velocity dispersion a via 



(1) 
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where 77 ps 10. The velocity dispersion a can be determined from high 
resolution spectroscopy and the half-light radius r^i from high resolution 
imaging. Photometry will yield the total luminosity L and the mass-to- 
light ratio M/L will provide an estimate of the underlying IMF. 

Using this technique Mengel et al. (2002) found that the IMF in 
selected SSCs in the Antennae galaxies is approximately consistent with 
Salpeter for a mass range of 0.1 — IOOMq. On the other hand, Smith & 
Gallagher (2001) found that in M82-F a lower mass cut-off of 2 - 3M 
is required for a MF with a Salpeter slope, i.e., M/L is 5x lower than 
normal. 

3.2 The problematic assumptions 

The approach described above depends on the validity of several as- 
sumptions: 

• The correct half-light radius rhi- Even gravitationally-bound, young 
clusters often deviate from radially symmetric, King-like cluster profile. 
A perfect example is R136, which has several bright cluster members at 
larger distance from the center. While those stars can be easily excluded 
from the light profile at the distance of the LMC, they blend with the core 
at larger distances (lower spatial sampling) and increase the measured 
rhi- This effect is illustrated in Figure 3. 

• A gravitationally bound cluster. While the two most massive HII re- 
gions in the Local Group, 30 Doradus in the LMC and NGC 604 in M 33, 
have about the same luminosity they have completely different morpho- 
logical structures. The core of 30 Dor, R136, is a massive compact cluster 
and presumably gravitationally bound, while NGC 604 consists of a sev- 
eral smaller, less luminous HII regions, which are not gravitationally 
bound together. However, at the distance of e.g., the Antennae galax- 
ies, where the scale of one WFPC2 pixel corresponds to about 9 parsecs, 
30 Dor and NGC 604 would be indistinguishable. 

• A constant M/L. Massive stars appear to be more concentrated toward 
the cluster center, see e.g., figure 16 in Brandl et al. (1996). We cannot 
distinguish here whether this is related to the star formation process 
(i.e., massive stars form in the densest regions), or the result of mass 
segregation. At any rate, the most massive stars dominate the light 
profile, and their concentration around the cluster center produces a 
gradient in M/L. 

While the above method certainly allows for rough estimates of dy- 
namical cluster masses, the combined uncertainties make it very difficult 
to provide sufficient evidence for anomalous IMFs in galaxies beyond the 
Local Group. 
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4. 



Method III: mid-IR fine-structure line ratios 



Luminous starbursts are often significantly extincted and require ob- 
servations at longer wavelengths. The mid-IR spectral range covered 
by ISO-SWS and Spitzer-IRS provides a wealth of spectral diagnostic 
features, although the spectra contain only spatially integrated informa- 
tion. 
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Figure 4- Model ratios of the mid- 
IR [Ne III] /[Nell] finestructure lines as 
a function of time and upper mass 
cut-offs. Figure from Rigby & Rieke 
(2004). 



Figure 5. Observed [Ne III] /[Nell] 
ratios from ISO-SWS for a variety 
of Galactic and extragalactic sources. 
The majority of luminous starburst 
galaxies fall below a ratio of one. Fig- 
ure from Thornley et al. (2000). 



4.1 The basic principle 

One of the most common diagnostics are the forbidden fine structure 
lines [NeIII]15.56/xm and [NeII]12.81/im. With an ionization potential 
of 63.45 eV, [Nelll] indicates the presence of early-type O-stars, while 
[Nell] with only 40.96 eV can be excited by less massive OB stars. Thus, 
the ratio of [Ne III] /[Nell] is a very sensitive measure of the strength of 
the radiation field, which depends on the presence of the most massive 
stars and on the age of the starburst, as illustrated in Figure 4. 

4.2 Problem: the model-observation discrepancy 

The observed line ratios (Figure 5) differ significantly from the val- 
ues predicted by the models such as STARBURST99 / CLOUDY (Lei- 
therer et al. 1999). While there seems to be a better agreement for 
low-met allicity galaxies and dwarf systems, the observed ratios of more 
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luminous starburst galaxies are often more than an order of magnitude 
lower than the theoretical predictions. 

Among the possible explanations being discussed is strong dependency 
on metallicity, a general lack of the most massive stars in intense bursts, 
or a larger age spread with complicated morphology of the starburst re- 
gion. Higher extinction can be ruled out as it would even increase the 
discrepancy (Figure 5). Rigby & Rieke (2004) proposed that very mas- 
sive stars in dense and metal-rich starbursts spend a significant fraction 
of their lifetime embedded in ultra-compact HII regions that prevent 
the formation of classical HII regions. Hence, the finestructure lines ex- 
pected from those embedded massive stars cannot even form, and the 
spectra would be dominated by the lower excitation [Ne II] line produced 
by the surrounding, less massive stars. 

In summary, mid-IR finestructure line diagnostics are an elegant and 
intriguing method to infer the relative amount of massive stars in a star- 
burst. However, given the complexity of a real starburst region and the 
model uncertainties claims of anomalous IMFs, just based on line ratios, 
need to be taken with care. This is clearly an area where significant 
improvements can be expected from Spitzer-IRS observations. 
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